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Abstract 

This study focuses on the multiple parton scatter- 
ing theory as a background for the new measurement 
of the inclusive cross section of the vector boson pair 
production at the LHC energies. The process under 
study is the double Drell-Yan annihilation. In this case, 
two quark- ant iquark annihilations occur independently 
in one proton-proton annihilation. The final state with 
two pairs of leptons (electron or muon pair) is investi- 
gated while the intermedial vector boson can be both 
gamma and Z. 

1. Introduction 

The traditional hadron collider experiment studies are 
usually based on the concept where one parton from ev- 
ery hadron takes place in the head-on collision of the 
colliding particles. With larger energies, the accelerated 
hadron is a composite object consisting of many partons 
from the sea. This increases the probability that two 
or more independent hard parton-parton scatterings can 
simultaneously appear in the same hadron-hadron in- 
teraction. This mechanism is generally called multiple 
parton scattering (MPS). 

In fact, there are two kinds of processes called multiple 
parton scattering, see Fig. [L]The disconnected scatter- 
ing was described in the previous paragraphs. During 
rescattering, one parton from the first-beam hadron in- 
teracts with two or more partons from the second-beam 
hadron. This process is significant for heavy ion colli- 
sions but is highly suppressed at proton-proton collisions 
and will not be further involved. 





Disconnected scattering 



Rescattering 



Fig. 1. Two kinds of multiple parton scattering: the discon- 
nected scattering and the rescattering. 



If all the parton interactions occur independently then 
the number of hard scatterings h at a given impact pa- 
rameter b and CMS energy ^/s of the given hadron- 
hadron interaction follows the Poissonian statistics: 



Ph(b,s) 



< n(b,s) > h o<n(M> 



hi 



(1) 



where < n(b,s) > is the average number hard scatterings 
and can be later expressed as (7). 

The additional parton scatterings proceed mostly via 
soft or semi-hard interactions at energies that have been 
reached up to now, but the chance of hard scattering will 
rise with the increase of the c.m.s. energy (like 7, 10, 
and 14 TeV planned for the LHC) and even faster than 
the classical single parton scattering production cross 
section. 



2. Double Parton Scattering 

The double parton scattering (DPS) is the simplest 
case of the MPS. Figure [271 shows two unspecified hard 
parton interactions taking part in one hadron-hadron in- 
teraction, and they are not connected or correlated in 
any kinematic observable. The processes are separated 
in transverse space by a distance of the magnitude of the 
hadron radius. The source of possible correlations would 
be an emission of virtual connecting (C) gluons by inter- 
acting partons. The C-gluons connect the hard processes 
(dotted lines in Fig. [2]), while the non-connecting (NC) 
gluons contribute to each hard process separately (curly 
lines in Fig. [2]). The CDF collaboration analyzed a 16 
pb 1 sample of pp — >• j/tt + 3 jets + X data pQ and did 
not find any apparent evidence of correlations in parton 
momentum fractions x, nor any evidence of kinematic 
correlations in mass, pr, or p z . 




Fig. 2. Schematic picture of double parton scattering with pos- 
sible virtual gluon connections. 



The multiple parton scattering processes offer a solu- 
tion for the unitarity violation of the total integrated 
cross section for hard collisions which is an already 
known fact appearing in the region near the perturba- 
tive threshold. The total cross section calculated from 
the parton model, assuming only single parton scatter- 
ing, predicts a constant rate in the whole perturbative 
transverse momentum range. However, there are many 
processes, especially multi-jet production, contributing 
over the predicted value, see Fig. [3] This increased ac- 
tivity can be explained by assuming the MPS. 



pTmin for pp at 14 TeV 




Fig. 3. The rise of the jet-production cross section over the total 
cross section at a low transverse momentum region. 



©2009 by Universal Academy Press, Inc. 



2 



For instance, DPS can produce the same final state 
(and with a comparable production rate) as the stan- 
dard single parton process with addition of perturbative 
corrections (multi-jets events) or by t-channel quark ex- 
change (IVB boson pair production). 

The importance of the MPS lies also in its ability to 
provide information on the spatial distribution of par- 
tons within the proton, as well as all other possible cor- 
relations among them. Its observation could bring a new 
parton concept complementary to the common descrip- 
tion, through the single-parton distribution functions. 

The DPS cross section can be factorized by the 
double parton distribution function T(xi,X2',b) as 
a function of two momentum fractions of the inter- 
acting partons and of their relative transverse distance b: 

— y J p c T r A (x 1 ,x 2 ',b)a(x 1 ,x' 1 )a(x 2 ,x 2 )r B (x / 1 ,x 2 ;b) 

x dx\dx' 1 dx2dx' 2 ^ (2) 

where m = 1 for indistinguishable partons processes and 
m = 2 for distinguishable partons processes. 

3. Hadron Matter Distribution 

The above mentioned DPS cross section formula [2] con- 
tains the double parton distribution functions. Current 
models, based on the inclusiveness of the standard par- 
ton distribution functions (PDF), allow the double par- 
ton distributions to be factorized by standard PDF's in 
convolution with the function of b expressing the spatial 
distribution of partons in the transverse plane and their 
overlap: 



T(x 1 ,x 2 ;b) = f(x 1 )f(x 2 )F(b). 



(3) 



The next step is to express the total cross section, 
where the number of parton interactions would appear 
as the index of the sum from one to infinity. Using 
the optical theorem and neglecting spin effects, the total 
cross section can be related to the imaginary part of the 
Fourier transform of the elastic scattering amplitude in 
the impact parameter a(b,s): 



°tot = 47r J d 2 blm (a(b, s)) . 



(4) 



The MPS model also stems from the eikonal approx- 
imation framework. Here, the elastic amplitude can be 
written down in terms of the eikonal function x(6, s) as 



a(b, s) 



-x(M) 



1 



2i 



(5) 



and similarly the inelastic part of the total cross section: 
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(6) 



The assumption of the independency of individual par- 
ton interactions leads to the Poissonian model, as men- 
tioned above. The average number of secondary hard 
scatterings is thus given by: 



< n(b,s) >= F(b)a hard (s). 



(?) 



The evaluation of the overlap function F(b) then de- 
pends on the model of the internal hadron structure. 
Figure explains the parameters used in the Eq. [8j 
where the overlap function is expressed as a convolution 
of the form factor distributions of two incoming hadrons: 



F(b) = J d 2 b'p A {b')p E 



(6-6'). 



(8) 




Fig. 4. Transverse overlap of the hadrons. 



The overlap function also satisfies the normalization con- 
dition: 



J 7rd 2 bF(b) = 1. 



(9) 



There are several possible models of hadron matter 
offering the expressions for the function p, represent- 
ing the parton matter distributions in transverse space. 
For example, the single Gaussian shape assumes the mo- 
mentum of individual partons to be distributed homoge- 
nously in the whole hadron, while the double channel 
Gaussian model predicts the hadron to contain a smaller 
hard core surrounded by a shell made up of softer par- 
tons. 

The question of the parton correlations inside the 
hadron is still unanswered and can include only the spa- 
tial constraints, as in the Poissonian model, or can con- 
tain all the possible correlations, such as in momentum, 
spin, and color of the diparton system, see e.g. [2]. 

4. Effective Cross Section 

The situation becomes more complex after involving 
the different behaviors of every kind of parton in the the- 
ory. The overlap functions need to carry two additional 
indices characterizing the two interacting partons. The 
formula for the DPS cross section (j2j) has to be rewritten 
into the form: 



ijkl 



a l i4 l Fi(b)F{(b)dH, 



(10) 



where as is the classical single parton scattering cross 
section. 

The integral over the impact parameter b is then 
marked as the geometrical coefficient O 



KJ ki 



= f Ft{b)Fi{b)d 2 b. 



(11) 



These geometrical coefficients have to be taken into ac- 
count when one expects different partons to have a dif- 
ferent distribution in transverse space inside the hadron. 

The simplest model ignores the differences in the ge- 
ometrical coefficients and approximates all the spatial 
correlations as well as the uncertainty in the transverse 
space parton distribution by the scale factor a e f f named 
effective cross section: 

°7f f = J [F(b)] 2 d 2 b. (12) 
Eventually the DPS cross section [2] can be simplified to: 



2 a, 



eff 



(13) 
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The prediction of the effective cross section value dif- 
fers according to the model of hadron structure. The 
first intuitive estimation uses the average hadron radius 
r w 0.6 fm. The a e ff would be approximately 30mb. 
The most promising model currently is the double Gaus- 
sian model that predicts a e ff « limb [31 . Table 1. sum- 
marizes the effective cross section measurements by ex- 
periments at CERN and Fermilab @], 0, ©, 0. 



Table 1. The effective cross section measurements. 







final state 


(T eff (mb) 


AFS, 1986 


63 


4jets 


« 5 


UA2, 1991 


630 


4jets 


> 8.3(95%C.L.) 


CDF, 1993 


1800 


4jets 


— iZ - i -5.4 


CDF, 1997 


1800 


7 + 3 jets 


= 14.5 ±1.71^ 



5. Double Drell-Yan Process 



longer tail. The average number of interactions has been 
found to be four. 

Other kinematical distributions seem to be similar to 
the single parton but are not investigated in detail in 
this stage of the study. There are few suspicions that 
the Monte Carlo implementation of the MPS model still 
needs to be tuned. The production cross section is not 
provided because the MPS is treated as the underlying 
event activity and can not be set as part of the produc- 
tion cross section calculation. 




Fig. 6. Distribution of total number of parton interactions for 
5000 proton-proton events at y/s = lOTeV. There are always 
at least two qq — > 7* jZ — > e + e~ interactions. 




Fig. 5. The DPS mechanism of 7*/^ pair production. 



The possibility of creation of a neutral vector boson 
pair via the double parton scattering mechanism has al- 
ready been suggested many years ago [8], [9]. In this 
process, the two independent s-channel quark- antiquark 
annihilations occur in one hadron-hadron interaction, see 
Fig. El 

On the basis of the previous theoretical description, 

the double Drell-Yan production cross section is pre- 

2 

dieted to be approximately <jddy — 2^77 ~ (J S-10~ 7 ~ 
0.1 fb, where only the electron decay channel is taken 
into account. The single parton production cross section 
was taken to be as ~ 1.1 rib. 

The motivation for this study was the promise of the 
cleanest signal in the detector [10] in comparison to the 
multi-jet processes, even at the cost of a small cross sec- 
tion. 

The source of the data for this preliminary analysis 
was the Herwig++ 2.3.0 Monte Carlo event generator 
[TT] . The model of MPS is already implemented in Her- 
wig++. The studied process consisted of the primary 
hard subprocess of type MEqq2gZ2ll and of the addi- 
tional hard subprocess of the same type with the multi- 
plicity equal to one. The full matrix element was used. 
The final state was created by the two electron-positron 
pairs. The MPS model of Herwig++ performs the calcu- 
lation of the number of sub-processes independently and 
adds the general QCD2 — )> 2 matrix element to them. 

The simulation was performed for proton-proton col- 
lisions at y/s = lOTeV. The distribution of the total 
number of parton interactions for 5000 proton-proton 
events is shown in Fig. [6] The shape of the distribution 
corresponds to the Poissonian distribution with a little 



6. Summary and Conclusions 

Multiple parton scattering was experimentally mea- 
sured by several collaborations and its existence was 
proven. The previous experiments aimed their attention 
at the multi-jet final states or at photon-plus-jets final 
state, where the production cross section is sufficiently 
large and the MPS processes contribute significantly es- 
pecially at a low pt region. The motivation to search for 
the other processes, like vector boson pair production, is 
increased by the approaching LHC era. 

The theoretical model describing MPS is basically de- 
veloped, but it suffers from many sources of uncertainty. 
Especially the internal parton correlations, which is one 
of the unknown effects that needs to be investigated fur- 
ther. MPS is also a unique tool for geometrical coefficient 
measurement. It may resolve the size of hard/soft com- 
ponent of the proton, and may bring about new many- 
body parton distribution functions. 

The preliminary predictions for the double Drell-Yan 
process shows a very low cross section for the lepton 
decay channel and it will not be measurable at LHC. For 
the same reason, it will not be a source of background 
for Higgs or for di-boson production. 
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